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Abstract 



We demonstrate, for the first time, successful S-matrix to potential inversion 
for spin one projectiles with non-diagonal Sj,> yielding a Tr interaction. The 
method is a generalization of the iterative-perturbative, IP, method. We 
present a test case indicating the degree of uniqueness of the potential. The 
method is adapted, using established procedures, into direct observable to 
potential inversion, fitting cr, iTn, T20, T%\ and T22 for d + alpha scattering 
over a range of energies near 10 MeV. The Tr interaction which we find is 
very different from that proposed elsewhere, both real and imaginary parts 
being very different for odd and even parity channels. 

PACS numbers: 21.30.-x, 13.75.Cs, 25.10,+s 



1 



Typeset using REVTeX 



The inverse scattering problem, i.e. the problem of determining a potential from the 
scattering matrix, has not until now been solved for the scattering of spin one projectiles. 
Spin one inversion is substantially different from that for spin zero or spin half because, even 
for a spin-zero target, the elastic scattering '^-matrix' is no longer a complex number but 
a matrix. When spin-1 projectiles scatter from a spin-zero target, the conserved quantities 
are total angular momentum j and parity 7T. When tt = (— two values of orbital 
angular momentum, I and I', are, in general, coupled by a tensor interaction. Calculating 
the S'-matrix S 3 U , therefore requires a coupled channel calculation and inversion from S 3 W to a 
potential containing a tensor term requires a coupled channel inversion. Here we demonstrate 
a spin-one inverse scattering procedure. It is incorporated into a direct data-to-potential 
inversion procedure which is used to analyse 2 H - 4 He scattering data at around 10 MeV. 
We find evidence for a substantial parity dependent tensor interaction. 

Coupled channel inversion is here carried out using a generalization of the iterative 
perturbative (IP) inversion procedure JJ-^j, which is successful for Sij — > V(r) + 1 • sV so (r) 
inversion for spin half projectiles. Features of the IP method (e.g. the ability to handle a 
range of energies simultaneously and to include Majorana terms for all potential components) 
all apply. 

We seek a potential with many components, each component being labelled with index 
p identifying central, spin-orbit or tensor terms, each real or imaginary. The number of com- 
ponents further doubles where, as essential for light nuclei, parity dependence is permitted. 
We assume that the tensor force is of the form V® = T R V®(r) = ((s • f) 2 - 2/3)V {t) (r) [§, 
where V® is assumed to have both Wigner and Majorana terms. 

The IP method commences from a 'starting reference potential', SRP and the potential is 
corrected iteratively. At each iteration, all components, indexed by p, are modified by adding 
a superposition a n v n{ r ) where functions v n (r) are members of an 'inversion basis' (which, 
if required, can be chosen differently for different p.) The amplitudes are determined 
at each iteration from linear equations, based on an SVD algorithm, which successively 
reduce a 2 = Y \£>t ~ ^k\ 2 - F° r eacn partial wave k, Sj, is the 'target' S'-matrix and S| is 
for the potential at the current iteration. Here label k is a single index which identifies 
non-diagonal elements of S 3 W as well as the energy when S J w (Ei) for energies Ei are 
simultaneously inverted. The physical basis of the equations leading to 0$ is the the linear 
response of S& to small perturbations in the potential This result applies to the non- 

diagonal as well as diagonal terms. For any given set of conserved quantum numbers, certain 
channels will be coupled by the nucleus-nucleus interaction and we use labels k, A,//, v for 
these channels. Thus the matrix element of the nucleus-nucleus interaction V between the 
wavefunctions for channels k and A, corresponding to integrating over all coordinates but r, 
will be written V K \(r). The increment AS K \ in the non-diagonal S-matrix which is due to a 
small perturbation AV K \(r), is 

AS ^ = §y. E % K {r)AV, u (r)Tp vX dr (1) 

where ip UK is the z/th channel (first index) component of that coupled channel solution for 
the unperturbed non-diagonal potential for which there is in-going flux in channel k (second 
index) only. The normalisation is ip UK — > S Kl/ Ii K — S UK Oi v where J; and 0\ are incoming 
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and outgoing Coulomb wavefunctions for orbital angular momentum /; there is no complex 
conjugation in the integral. Starting from Eq.|l], spin-one inversion becomes a straightforward 
generalisation of the procedure described in Refs. || and is now implemented in the code 
IMAGO [|. 

Testing spin-1 S —>■ V inversion. We established that in forward (rather than inverse) 
mode, a given Tr potential led to the same S 3 U , and observables (including tensor analysing 
powers) as the standard deuteron scattering code DDTP ||. A convenient property of 
the IP inversion method is the fact that when a sequence of iterations converges to small 
a 2 , as defined above, we can be sure that we have found a potential which reproduces the 
set of target S 3 W to a precision which corresponds to that value of a 2 . It follows that the 
forward-mode test verifies that we have found a correctly defined potential. 

We have also tested the prescription for spin-1 S — > V inversion using a known potential 
to define the target S{ v . One set of results, for 20.2 MeV deuterons scattering from 4 He is 
shown in Fig. 1. The target potential |7J (solid line) has a complex central term and real 
1 • s and V® terms but no Majorana components. This potential reproduces the tensor 
analysing powers but with a Tr interaction much stronger than folding models || predict. 
The two inversion solutions shown in Fig. 1 both arise from a starting potential (SRP) with 
only one non-zero component, a real central Woods-Saxon potential 60 MeV deep. Although 
few partial waves contribute significantly to the inversion, only 3-4 iterations were necessary 
to give solution A (dotted) and solution B (dot-dashed), which reproduced all four volume 
integrals to better than 1 %, required not many more. Defects in the tensor and spin-orbit 
terms for r ~ reflect the fact that both 1 • s and Tr have zero diagonal matrix elements for 
I = 0. A test case involving parity dependence would necessitate a smaller inversion basis, 
and hence reproduce each component with less accuracy, since S 3 U , for the same range of j 
must determine twice as many components; this is relevant to the case described below. 

Direct data-to-potential inversion. The preferred method of applying inversion 
methods to experimental scattering data is to convolute the phase shift fitting with the 
S — > V inversion to form a direct data-to-potential inversion procedure. This can be coded 
in a single data fitting program. The two step procedure in which one first fits the data with 
S 3 U , and then subsequently inverts S 3 U , does not guarantee that the S 3 W for different energies 
are smoothly related like S 3 U , derived from a single potential. The direct inversion formalism 
has been described and applied to cases where there is no channel coupling and it 

carries over to the present case with no change in basic principle. It has been implemented 
in code IMAGO and the convolution was thoroughly tested through explicit checking of 
all the derivative terms which are required [PjLOfl. 

Application to deuteron scattering. We have analysed the scattering of tensor 
polarised deuterons from 4 He over a series of energies near 10 MeV using direct data-to- 
potential inversion. This case is of particular interest both because it should become fea- 
sible to do quite realistic calculations using RGM or similar methods and also because the 
deuteron is the archetypal easily polarisable halo nucleus. From the data set for d - alpha 
scattering tabulated in Ref. [ I I], we have selected data of Jenny et al \T2\ which is of high 
quality and covers a wide angular range. We fit the following observables, a, iTn, T 20 , 
T 2 i and T 22 , for the energies, 8, 9, 10, 11, 12 and 13 MeV. This range includes broad 2 + 
resonances and a region of strong mixing between the 1 + channels. Ref [13| applies direct 
data-to-potential IP inversion to data from the same source but without fitting T 2i . 
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We exploit a property of the IP method to fit all of the above data, for all the stated en- 
ergies, with a single potential. The potential is of fixed radial form but the overall strength of 
each of the imaginary components is proportional to the energy above the inelastic thresh- 
old. To minimize the number of parameters the real components are taken to be energy 
independent. 

It is well known that fits of ^-matrix to data are inevitably subject to discrete and 
continuous ambiguities and direct data-to-V inversion is not exempt. The uncertainties 
arising at the S — > V stage are less important. Various strategies are possible, such as the 
inclusion of a priori information from theory [[Hfl . Here we adopt the approach of using 
extremely small inversion bases (two or three Gaussian functions) so that, in effect, we find 
the smoothest, minimally energy-dependent, potential fitting the data at all six energies. 
We find that we can allow the spin orbit terms to be real and parity independent but the 
central and Tr tensor terms must be complex and and parity dependent. Inversion studies 
involving fits to both empirical data and RGM theory (T3.T^| show that interactions between 



light nuclei are parity dependent for energies up to tens of MeV per nucleon. As in earlier 
work, we present the real central term in the form of Wigner and Majorana components: 
Vw + (— 1)'Vm- However, V"w and the imaginary terms were treated somewhat differently: 
the data was reproduced most efficiently with inversion bases covering different radial ranges 
for the even and odd parity components. In Fig. 3 below we shall see that the odd parity 
y{t) terms cover a wider radial range. This is consistent with theoretical arguments [|17] for 



the existence of tensor interactions of quite different strengths and shapes for odd and even 
parity. The SRP consisted only of real, Wigner, central and spin-orbit terms. Convergence 
was very rapid and typically about 3 iterations were required. 

Figs 2 and 3 show two alternative potentials fitting the data for the six energies. The 
energy dependent imaginary parts were evaluated at 10 MeV and Fig. 4 shows the quality 
of fit to the 10 MeV data. The fit to the data at the ends of the energy range was poorer 
suggesting that a somewhat less simple energy dependence is required. The y^/N values of 
8.59 and 8.08 are for the data at all six energies.^ 

Different components of the potential are determined to different degrees of uniqueness. 
The real central and spin-orbit Wigner potentials shown in Fig. 2 are well determined and 
consistent with predictions of RGM theory. The imaginary central component is highly 
parity dependent and is emissive near the origin. All potentials giving a good fit to the 
data over this energy range have such an emissive region near the nuclear centre. Such 
emissive regions are commonly present [Hjj in potentials found by inverting RGM S'-matrix 



elements which manifestly do not break unitarity. Emissive regions often appear in local 
potentials which represent exchange and channel coupling non-locality. In the present case, 
the unitarity limit is only broken, to a small degree, for L = 2, J = 2 for the y 2 /N = 8.59 
and is not broken for the better fit. 

The V®(r) interaction shown in Fig. 3 departs markedly from previous phenomenology. 
The even parity real and imaginary terms are very strong near r = 1. This is significant for 



1 For simplicity in this initial study, we made no selection, filtration or normalisation of input data 
as is usual in phase shift analyses. 
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r > about 1.0 fm; for r < 1 the even-parity V®(r) interaction is ill-determined since the 
diagonal matrix element of Tr is zero for I < 2 so the comment concerning irregularities at 
the nuclear centre in Fig. 1 applies here. The combination of strong, short-ranged even- 
parity components and weaker, long ranged odd-parity components is precisely the pattern 
predicted to arise from a deuteron exchange mechanism discussed in Ref. |17| ; see also 



Refs. [fl8| , |T9"|| which also present a V"^ interaction strongly peaked inside r = 2 fm, but of 
the opposite sign. 

The present analysis firmly establishes that most terms in the d - 4 He potential are 
substantially parity dependent, casting doubt on parity independent d - 4 He potentials, 
including the potential of Frick et al with its remarkably large V^' interaction. RGM 
theory [0,|1^] predicts parity dependence for central components, but the strong parity 
dependence in V"w suggests a new physical process. 

In conclusion: we have demonstrated for the first time S — > V inversion for spin-1 
projectiles leading to a tensor interaction which couples channels. This has been incorporated 
into a direct data-to-potential inversion procedure yielding a remarkable new kind of Tr 
interaction. The general procedure could be applied to other cases with channel spin one 
such as p - 3 H scattering. Spin-1 S — > V inversion makes it possible to study theories for 
dynamic polarization and exchange contributions to central and non-central forces for spin-1 
projectiles. The results can then be directly related to phenomenology since, as we have 
shown here, the inversion procedure can be made part of a powerful, rapidly converging, 
phenomenological method which, at low computational cost, fits multi-energy datasets with 
a single, energy dependent potential. 
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FIGURES 



FIG. 1. For deuterons scattering from 4 He at 20.2 MeV, the full line represents the potential 
of Frick et al, the long dashes represent the SRP (non-zero only for the real central term), the dots 
represent inversion solution A, and the dot-dashes represent solution B. 

FIG. 2. From top: Real, Wigner and Majorana, central components; real, Wigner spin-orbit 
component; imaginary, central, even and odd parity components. All were evaluated at 10 MeV, 
and correspond to inverted potentials with % 2 /N = 8.59 (solid line) and x 2 /N = 8.01 (dashed 
line) . 

FIG. 3. Tensor components of the potentials of Figure 2, from top: Even and odd, real then 
even and odd imaginary. Again, solid line is for x 2 /N = 8-59 and dashed line is for x 2 /iV = 8.01. 

FIG. 4. For deuterons scattering from 4 He at a laboratory energy of 10 MeV, u(9), iTn(9), 
T2q{9), T2\{9), T22{9) given by the potentials shown in Figures 2 and 3, with x 2 /-/V = 8.59 (solid 
line) and x 2 /N = 8.01 (dashed line), compared with the data of Jenny et al (solid points). 
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Fig. 3 
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